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ABSTRACT 

Using  a  model  for  the  heat  flux  through  a  single 
drop,  a  computer  simulation  of  the  dropwise  condensation 
process  was  carried  out.   For  a  given  set  of  condensing 
conditions,  the  computer  model  predicts  heat  transfer 
coefficients  and  drop  distributions  which  agree  with 
experimental  observations.   Subsequently,  an  analysis  was 
conducted  with  the  model  by  varying  the  condensing  para- 
meters of  site  density,  saturation  temperature,  and  tempera- 
ture difference.   From  a  comparison  of  the  analysis  with 
experimental  results,  the  effects  of  the  previously 
mentioned  parameters  were  determined. 
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NOMENCLATURE 

a  =  constant  defined  in  Eq.  (A. 15) 

A  =  Area  Cft  ) 

b  =  constant  defined  in  Eq.  (A. 15) 

c  =  constant  defined  in  Eq.  (A. 15) 

d  =  constant  defined  in  Eq.  (A. 15) 

D  =  drop  diameter  (ft) 

h   =  heat  transfer  coefficient  defined  in  Eq.  (F.3) 

CBtu/hr  ft2  °f) 
H   =  latent  heat  of  vaporization  (Btu/lbm) 

K  -  thermal  conductivity  (Btu/hr  ft°F) 

M  =  molecular  weight 

Q  =  rate  of  heat  transfer  (Btu/hr) 

r  =  radius  of  drop  (ft) 

R  =  resistance  (hr  °F/Btu) 

R  -   universal  gas  constant 

t  =  time  or  drop  age  (seconds) 

T   =  temperature 

3 

v   =  specific  volume  (ft  /lbm) 

3 

V   =  volume  of  drop  (ft  ) 

a   =  mass  transfer  accommodation  coefficient, 
condensation  coefficient  (dimensionless) 

AT  =  temperature  difference  (°F) 

3 
p   =  density  (lbm/ft  ) 

a   =  surface  tension  (lbf/ft) 
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NOMENCLATURE 
(Cont'd) 


SUBSCRIPTS 
c   =  curvature 
cyl  =  cylinder 
dc  sa  drop  conduction 

fg  =  refers  to  change  by  condensation 
g   =  vapor 
i   =  interfacial 
max  =  maximum 


mm  = 

minimum 

nuc  = 

nucleation 

s   = 

saturation 

t 

total 
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CHAPTER  I 
INTRODUCTION 
Over  the  past  few  decades  condensers  have  remained 

essentially  unchanged,  while  other  components  of  the  steam 
power  plant  have  been  improved.   The  potential  for  improvement 
of  the  condenser  exists  if  the  phenomenon  of  dropwise  con- 
densation can  be  utilized.   Dropwise  condensation  is  the 
process  where  vapor,  condensing  on  a  cold  and  non-wetting 
surface,  forms  liquid  drops  which  grow  and  coalesce  and 
ultimately  roll  off.   For  a  surface  which  the  liquid  wets, 
a  continuous  liquid  film  forms  and  covers  the  surface. 
This  latter  process  is  called  filmwise  condensation  and  is 
the  basis  of  current  condenser  design.   Experimentally, 
heat  transfer  coefficients  for  dropwise  condensation  have 
been  found  to  be  an  order  of  magnitude  larger  than  the 
coefficients  for  filmwise  condensation  [1] .   With  this 
superior  heat  transfer  property,  condensers  can  be  made 
smaller  by  weight  and  volume. 

There  are  three  problems  retarding  the  use  of 
dropwise  condensation.   First,  there  is  today  no  well  proven 
surface  or  surface  coating  which  promotes  the  formation  of 
drops.   Secondly,  the  presence  of  noncondensable  gases  in 
the  vapor  reduces  the  heat  transfer  coefficient.   There  is  a 
need  for  a  means  of  reducing  the  noncondensables  to  an 
acceptably  low  limit.   Finally,  there  is  a  lack  of  understand- 
ing of  the  effect  of  certain  condensing  parameters.   This 
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study  is  to  aid  understanding  by  providing  a  theoretical 
prediction  of  the  heat  transfer  properties  of  dropwise 
condensation  by  use  of  computer  simulation  of  the  growth  and 
coalescence  of  drops  from  nucleation  to  the  size  at  which 
the  drops  slide  off  the  surface. 

This  investigation  models  condensation  on  a 
vertical  surface.   Initially  the  surface  is  bare.   Then 
steam,  without  noncondensable  gases,  is  exposed  to  the 
surface.   Drops  form  at  nucleation  sites  which  are  probably 
randomly  distributed.   They  are  estimated  to  be  on  the  order 
of  one  hundred  million  sites  per  square  centimeter  [2] .   The 
nucleation  size  drops  are  of  the  order  of  one-tenth  of  a 
micron  in  diameter.   These  small  drops  initially  are  distri- 
buted over  the  entire  surface  and  begin  growing  very 
rapidly  due  to  condensation  of  vapor.   When  neighboring  drops 
touch,  they  coalesce  into  a  single  new  drop.   If  a  nucleation 
site  is  exposed  by  the  coalescence  activity,  a  new  nucleation 
size  drop  appears.   The  process  of  growth  and  coalescence 
continues  until  a  drop  reaches  a  size  at  which  the  surface 
tension  forces  holding  the  drop  on  the  surface  are  exceeded 
by  the  gravity  force  on  the  drop  mass.   Drops  of  the  depart- 
ing size  are  observed  to  slide  down  the  surface  and  sweep 
other  drops  with  them.   This  departing  action  also  exposes 
nucleation  sites. 
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In  1930  Schmidt,  Schurig,  and  Sellschopp  [3]  first 
investigated  dropwise  condensation  as  a  means  of  heat  transfer 
Since  their  investigation,  there  has  been  a  controversy  as  to 
the  origin  of  the  drops  and  as  to  whether  the  heat  is 
transferred  through  the  drops  or  through  the  bare  surface 
between  the  drops .  The  opinions  can  be  divided  into  two  broad 
theories.   First,  the  "film  theory"  contends  that  vapor 
condenses  on  the  bare  surface  between  the  drops  and  that  the 
resulting  liquid  forms  into  drops.   Variations  of  this  theory 
are  supported  by  the  work  of  early  investigators  such  as 
Jacob  14] ,  Eucken  [5] ,  and  Emmons  [6] .   More  recent  investi- 
gators have  proposed  the  "nucleation  theory"  which  disagrees 
with  the  "film  theory."   The  proponents  of  the  "nucleation 
theory"  contend  that  drops  nucleate  from  discrete  and 
randomly  distributed  sites  on  the  condensing  surface. 
Furthermore,  the  theory  states  that  there  is  no  liquid 
between  the  drops  and  that  the  path  of  heat  transfer  is 
through  the  drops  only.   In  1935  Tammann  and  Boehme  [7] 
were  the  first  of  several  to  observe  that  drops  form  at  the 
same  sites  on  successive  condensation  cycles.   In  1964  Umur 
and  Griffith  [8]  proved  that  there  is  no  liquid  between  the 
drops.   Using  thermodynamic  considerations  and  an  optical 
study,  they  demonstrated  that  a  film  no  greater  than  a 
monolayer  thick  could  exist  between  the  drops  on  the  surface. 
Umur  and  Griffith's  work  is  final  strong  proof  of  the 
validity  of  the  "nucleation  theory." 
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In  agreement  with  the  nucleation  theory,  Rose 
[9,10]  was  the  first  to  develop  a  method  of  predicting  the 
heat  transfer  pei formance  of  dropwise  condensation  by 
analyzing  the  heat  transferred  through  a  single  drop  and  then 
calculating  the  heat  flux  through  the  condensing  surface  by 
assuming  a  drop  distribution.   By  fitting  four  constants  into 
his  equations,  Rose  was  able  to  predict  values  which  were 
comparable  to  experimental  results.   Continuing  with  a 
similar  approach,  Graham  [2] ,  further  developed  an  analytical 
expression  for  the  heat  flux  through  a  drop  which  included 
the  heat  transfer  resistances  of  drop  curvature,  interfacial 
resistance  between  the  liquid  and  vapor,  and  conduction.   In 
order  to  authenticate  his  model,  he  experimentally  measured 
the  steady  state  drop  distribution.   The  smallest  drops 
which  Graham  measured  and  counted  were  ten  microns  in 
diameter.   Drops  of  smaller  diameter  could  not  be  accurately 
observed  because  of  a  resolution  limitation  of  the  optics. 
The  small  drops  are  moving  very  fast  and  have  diameters  on  the 
order  of  the  wave  length  of  visible  light.   By  combining  his 
experimental  distribution  and  his  analytical  model,  Graham 
found  that  fifty  percent  or  more  of  the  heat  flux  was 
transferred  through  drops  smaller  than  ten  microns  in 
diameter  and  concluded  that  his  heat  transfer  model  could 
not  be  verified  unless  the  drop  distribution  was  known 
over  the  entire  range  of  drop  sizes.   Since  a  complete  drop 
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distribution  cannot  be  experimentally  determined,  Graham 
recommended  computer  simulation  of  the  dropwise  condensation 
process  as  a  means  of  obtaining  the  complete  drop  distribution. 

The  first  computer  simulation  of  steady  state 
dropwise  condensation  was  carried  out  by  Gose,  Mucciardi, 
and  Baer  [11] .   The  model  included  nucleation,  drop  growth  due 
to  condensation,  coalescence,  and  removal  of  drops  of  the 
departing  size.   Their  expression  for  drop  growth  from  con- 
densation considered  only  the  conduction  resistance  to  the 
heat  flux  through  the  drop.   For  a  departing  drop  with  a 
diameter  of  two  millimeters  and  for  a  site  density  of 
5  x  10   sites  per  square  centimeter,  their  model  predicted 
a  heat  transfer  coefficient  of  approximately  460  British 
Thermal  Units  per  foot  square  hour  degree  Fahrenheit.*  This 
value  is  more  than  an  order  of  magnitude  lower  than  experi- 
mentally obtained  coefficients. 

A  later  study  by  Tanasawa  and  Tachibana  [12]  used 

a  similar  technique  and  predicted  a  coefficient  of  3,700  for 

3 
a  site  density  of  4  x  10  .   This  prediction  is  lower  than 

the  value  of  40,000  measured  by  Graham  [2]  at  the  same 

condensing  conditions.   Both  of  the  previous  computer 

simulations  modeled  drops  of  the  entire  range  of  sizes  on 

a  fixed  area.   This  approach  and  their  computer  facilities 

*For  this  thesis  heat  transfer  coefficients  will  have  the 
dimensions  of  British  Thermal  Units  per  foot  square  hour 
degree  Fahrenheit  and  site  densities  will  be  expressed  as 
sites  per  square  centimeter. 
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limited  their  investigations  to  artifically  low  site 
densities.   Gose  suggested  that  theoretical  values  would  more 

closely  match  experimental  results  if  larger  site  densities 

7       8 
of  10   to  10   could  be  modeled.   These  two  computer  simula- 
tions and  Graham's  work  provided  the  essential  background 
for  this  thesis. 

Using  Graham's  model  of  heat  flux  through  a  single 
drop,  a  computer  simulation  of  the  dropwise  condensation 
process  was  carried  out  for  this  thesis  in  three  parts. 
First,  a  computer  model  was  constructed  which  predicts  heat 
transfer  coefficients  and  drop  distributions  for  a  given  set 
of  condensing  conditions.   Secondly,  a  parametric  analysis 
was  conducted  varying  the  significant  parameters,  site 
density,  saturation  temperature  of  the  vapor,  and  the  tempera- 
ture difference  between  the  vapor  and  the  condensing  surface. 
The  third  and  final  part  of  this  thesis  consisted  of  a  compari- 
son of  experimental  results  and  model  predictions. 
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CHAPTER  II 
COMPUTER  SIMULATION 

The  Computer  Model 

To  model  a  representative  condensing  surface  of  one 
square  centimeter  would  require  the  use  of  approximately  one 
hundred  million  nucleation  sites.   The  computer  storage 
necessary  to  identify  the  location  and  activity  of  these 
nucleation  sites  was  not  available  with  the  IBM  360-65  Model 
Computer  with  a  Fortran  G  Compiler  used  for  this  study.   If 
core  storage  were  available,  the  computational  time  for  high 
site  densities  would  be  large  and  the  procedure  would  be 
uneconomical.   To  solve  the  storage  constraint  and  to  build 
an  accurate  and  efficient  computer  simulation,  an  innovative 
modeling  technique  was  developed  to  model  the  dropwise  con- 
densation process  on  an  initially  drop  free  surface. 

The  simulation  models  a  single  cycle  in  which  a  drop 
of  the  departing  size  is  formed  on  an  initially  bare  surface. 
To  start  the  modeling  process,  the  surface  is  exposed  to  an 
unlimited  supply  of  vapor.   Small  primary  drops  appear  at 
randomly  distributed  nucleation  sites.   These  drops  increase 
their  size  by  condensation  of  vapor  on  their  surface  and  by 
coalescence  with  neighboring  drops.   As  the  clock  time 
increases,  the  average  drop  size  existing  on  the  surface 
increases.   Eventually,  a  drop  of  the  departing  size  is 
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produced  and  slides  off  the  surface.   The  simulation  process 
stops  when  the  first  drop  of  departing  size  is  formed. 

The  approach  of  the  simulation  is  one  of  "viewing "' 
the  condensing  surface  with  decreasing  magnification  as  the 
average  drop  size  on  the  surface  increases.   With  a  micro- 
scope one  can  see  smaller  objects  by  increasing  the 
magnification,  but  by  doing  so,  the  field  of  vision  is 
decreased.   Similarly,  the  first  area  of  interest  is 
"viewed"  with  the  highest  magnification  and  is, therefore,  the 
smallest  area  modeled.   As  the  clock  time  increases,  the 
average  drop  size  increases.   When  the  area  of  interest  is 
no  longer  representative,  the  magnification  is  decreased  and 
hence  a  larger  area  is  visible.   Ensuing  drop  growth 
necessitates  lower  and  lower  magnifications  until  departing 
drops  are  visible.   An  area  of  interest  for  a  given  magnifi- 
cation is  called  a  stage.   Table  1  shows  how  the  stages  are 
arranged  for  a  representative  calculation.   Drops  of  the 
smallest  size  are  not  visible  on  later  stages,  but  their 
presence  is  accounted  for. 

The  modeling  technique  starts  with  two  assumptions. 
The  simulation  models  the  continuous  process  with  discrete 
time  intervals  or  "time  steps".   Drop  growth  from  condensation 
and  drop  coalescence  are  modeled  to  occur  sequentially  within 
the  time  step.   Whenever  the  simulation  clock  time  is  increased 
by  the  time  step,  a  routine  to  determine  drop  growth  is 
carried  out,  then  the  coalescence  activity  is  modeled  by  a 
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separate  routine.   This  procedure  is  considered  to  give  an 
accurate  representation  of  the  actual  continuous  process. 
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TABLE  1 
ARRANGEMENT  OF  STAGES 


Stage 
Number 

Area 

in 

Square  Microns 

1 

103 

2 

104 

3 

105 

4 

106 

5 

107 

6 

3. 

24  x  108 

Range  of  Radii  of  "Visible" 
Drops  in  Microns 

0.10  to  3.75 

1.59  to  12.04 

8.03  to  40.64 

27.10  to  137.17 

91.45  to  462.96 

308.64  to  1250.00 


Note:   This  stage  arrangement  is  for  a  saturation  temperature 

g 

of  212°F  and  for  a  nucleation  site  density  of  10 

sites  per  square  centimeter. 
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The  growth  routine  determines  a  drop's  size  from 
an  expression  of  drop  age  as  a  function  of  drop  diameter. 
The  age  is  increased  by  the  time  step.   The  drop's  new  size 
due  to  condensation  during  the  time  step  is  calculated  using 
a  relationship  for  drop  diameter  as  a  function  of  drop  age. 
A  detailed  discussion  of  the  growth  routine  is  included  in 
Appendix  A. 

Upon  completion  of  the  growth  routine,  a  check  for 
coalescence  is  carried  out.   The  coalescence  routine  employed 
by  this  study  is  innovative  and  consists  of  a  method  of 
selectively  searching  for  coalescences  between  drops  and  for 
the  covering  and  uncovering  of  nucleation  sites.   Each  drop 
is  considered, in  turn, from  left  to  right  on  the  stage  area  as 
the  focal  point  of  the  search.   A  search  square  of  two  drop 
diameters  on  a  side  is  centered  on  the  drop  of  interest. 
Only  nucleation  sites  and  drops  with  centers  inside  the 
search  square  are  considered.   A  determination  is  made 
whether  or  not  the  central  drop  blankets  any  of  the  active 
sites  or  coalesces  with  any  of  its  neighbors.   If  an 
uncovered  site  is  found  to  be  blanketed,  it  is  made  inactive. 
If  the  drop  of  interest  is  found  to  be  touching  a  neighbor 
drop  of  the  same  size  or  smaller  inside  the  square,  they 
coalesce  to  form  a  single  new  drop  at  the  central  drop's 
location.   After  the  coalescence  occurs,  a  check  for  the 
uncovering  of  inactive  sites  is  made.   Sites  uncovered  by  the 
coalescence  are  made  active. 
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The  coalescence  has  increased  the  size  of  the 
central  drop.   Consequently,  an  enlarged  search  area  must  be 
considered  so  that  all  possible  coalescences  of  the  central 
drop  with  neighboring  drops  of  the  same  size  or  smaller  are 
found.   After  each  coalescence,  the  search  area  is  enlarged. 
When  the  central  drop  fails  to  coalesce,  the  search  recommen- 
ces around  the  next  drop  to  the  right  on  the  stage  area. 
When  all  drops  have  been  the  center  of  a  search  area,  the 
coalescence  routine  ends. 

The  central  drop  only  coalesces  with  drops  of  the 
same  size  or  smaller.   By  following  this  procedure  for  all 
search  squares,  all  possible  coalescences  are  found.   This 
search  technique  has  the  advantages  of  requiring  a  small 
amount  of  core  storage  and  is  reasonably  efficient  for  one 
thousand  or  fewer  sites  per  stage  area.   Appendix  B  contains 
a  more  detailed  explanation  of  the  coalescence  routine. 

To  start  the  first  stage  requires  selection  of  a 
time  step  value.   If  the  time  step  is  too  large,  the  drops 
nucleated  on  the  first  time  step  will  overlap  by  a  large 
amount.   Drops  are  observed  to  coalesce  upon  first  touching. 
Significant  overlapping  is  a  serious  departure  from  nature 
at  high  site  densities.   The  time  step  for  a  given  site 
density  should  be  small  enough  to  produce  primary  drops 
which  do  not  artifically  overlap. 
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The  time  step  can  also  be  too  small.   The  problem 
arises  because  the  approximation  for  drop  radius  as  a  function 
of  age  becomes  inaccurate  for  very  small  drop  radii.   For  this 
reason,  the  time  step  should  be  large  enough  to  produce 
primary  drops  which  are  three  to  four  times  the  minimum  drop 
radius  CAppendix  A) . 
The  First  Stage 

o 

For  a  site  density  of  10  ,  the  first  stage  has  one 
thousand  randomly  distributed  nucleation  sites  on  a  square  of 
thirty-one  microns  on  a  side.   On  the  first  time  step  of  the 
first  stage,  the  growth  routine  causes  hemispherical  drops  of 
uniform  size  to  appear  centered  on  nucleation  sites.   The  size 
of  the  uniform  drops  is  commensurate  with  the  time  step. 
Next,  the  coalescence  routine  is  carried  out.   If  neighboring 
drops  touch ,  they  coalesce  to  form  a  single  new  drop  centered 
on  the  nucleation  site  of  the  formerly  larger  drop.   At  the 
end  of  the  coalescence  routine,  the  drops  are  no  longer  of 
uniform  size.   Their  sizes  vary  depending  on  the  number  of 
coalescences  which  they  have  undergone.   The  stage  surface 
contains  nucleation  sites  which  were  uncovered  by  the 
coalescence  activity  as  well  as  drops  of  different  sizes. 
Other  nucleation  sites  are  inactive  because  they  are 
covered  by  a  drop.   Upon  completion  of  the  coalescence 
routine,  the  simulation  continues  by  increasing  the  clock 
time  by  the  time  step. 
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During  the  second  time  step,  new  drops  nucleate 
and  grow  on  the  nucleation  sites  uncovered  during  the 
previous  coalescence  routine.   Drops  already  existing  on  the 
surface  grow  by  condensation.   The  amount  of  growth  depends 
on  the  drop's  size  and  the  magnitude  of  the  time  step.   The 
coalescence  check  is  made,  and  the  clock  is  again  increased. 
The  growth  and  coalescence  sequence  with  increasing  time 
continues  until  the  stage  termination  criteria  is  reached. 

The  stage  termination  criteria  prevents  "border 
effects".   Drops  near  the  border  or  edge  of  the  stage  area 
have  an  abnormally  low  number  of  neighbors;  therefore,  they 
have  less  opportunity  to  increase  their  size  by  coalescence. 
The  effect  of  this  unrepresentative  activity  is  minimized  by 
modeling  drops  on  a  given  stage  only  until  they  become 
relatively  large  when  compared  to  the  surface  area.   There 
are  two  tests  in  the  program  for  stage  termination.   If  a 
single  drop  is  formed  which  covers  greater  than  ten  percent 
of  the  stage  area,  the  stage  ends.   A  category  of  smaller 
drops  may  also  cause  termination  if  they  cover  more  than 
twenty  percent  of  the  stage  area.   These  drops  have  radii 
which  are  forty  percent  smaller  than  the  radius  of  the  single 
drop  which  can  cause  termination. 

The  amount  of  heat  transferred  through  the 
condensing  surface  can  be  determined  from  the  amount  of 
liquid  contained  in  the  drops  on  the  surface.   The  amount  of 
heat  transferred  during  a  time  step  is  simply  the  total 
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volume  change  of  all  drops  on  the  surface  multiplied  by  the 
density  of  the  liquid  and  the  latent  heat  of  transformation. 
Since  the  rate  of  heat  transfer,  condensing  surface  area, 
and  temperature  difference  between  the  vapor  and  the  surface 
are  known;  the  heat  transfer  coefficient  for  each  time  step 
can  be  calculated.   From  the  total  amount  of  liquid  mass  on 
the  surface  at  a  given  clock  time,  the  time  averaged  heat 
transfer  coefficient  for  the  process  can  be  obtained  and 
is  the  primary  output  of  the  stage.   Figure  1  illustrates 
the  time  averaged  heat  transfer  coefficient  as  a  function  of 
time. 
The  Second  Stage 

Since  the  area  under  study  in  the  second  stage  is 
much  larger  than  the  first  stage,  the  nucleation  cannot  be 
modeled  in  detail  as  it  was  on  the  first  stage.   Rather,  it 
is  assumed  that  on  the  bare  area,  nucleation  and  growth  for 
the  time  interval  equal  to  t,  results  in  the  same  heat  flux 
per  unit  area  as  growth  of  the  first  stage  from  time  zero  to 
time  t,.   Thus  the  heat  flux  for  the  bare  area  is  calculated 
as  the  product  of  bare  area,  time  interval  t,  ,  and  average 
heat  transfer  coefficient  at  time  t,;  see  Figure  1.   This 
heat  flux  results  in  the  creation  of  an  amount  of  liquid  on 
the  bare  area.   Drops  of  radius  r2  are  placed  on  uncovered 
locations  on  the  bare  area  to  account  for  this  amount  of 
liquid  for  each  time  step  of  the  second  stage. 
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The  starting  conditions  of  the  second  stage  are 
determined  by  an  interative  matching  procedure  which  matches 
the  results  of  the  first  stage  with  the  start  up  of  the 
second  stage.   The  start  of  the  second  stage  is  at  the  time 
t,  Csee  Figure  1)  which  is  half  the  cycle  time  of  the  first 
stage.   The  initial  drops  of  the  second  stage  are  assumed  to 
be  uniform  in  size.   A  drop  radius,  r2,  is  assumed.   The 
number  of  drops,  N2 ,  is  calculated  so  that  the  total  amount 
of  liquid  on  the  surface  at  time  t,  is  the  same  for  the  first 
and  the  second  stage;  this  makes  the  average  heat  transfer 
coefficient  of  the  first  and  second  stage  equal. 

The  N2  drops  are  distributed  on  a  uniformly  spaced 
net  of  locations,  the  number  of  locations,  N,  ,  being  larger 
than  N2«   At  each  location  a  random  number  is  used  to  deter- 
mine if  a  drop  is  present.   All  N~  drops  are  located  on  the 
NL  locations.   The  N.  locations  are  uniformly  spaced  to 
ensure  that  none  of  the  drops  coalesce  immediately  after 
they  are  placed  on  the  surface. 

After  all  of  the  drops  have  been  placed  on  the 
surface,  the  time  is  increased  to  2t,  and  the  drops  grow  by 
condensation.   Additional  drops  of  radii  r2  are  placed  on 
uncovered  locations  to  account  for  the  condensing  productivity 
of  the  invisible  drops  on  the  apparently  bare  area.   The 
average  heat  transfer  coefficient  at  time  2t,  is  calculated 
and  compared  to  the  value  found  from  the  first  stage.   If  the 
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values  disagree,  a  new  radius  r~  is  assumed  and  the  process 
is  repeated  until  the  average  heat  transfer  coefficient  of 
the  first  and  second  stages  agree  at  times  t,  and  2t,. 

When  the  proper  value  of  r~  is  obtained,  the  second 
stage  is  allowed  to  proceed  with  time  steps  equal  to  t-,  .   The 
process  of  growth  and  coalescence  of  existing  drops  and  the 
creation  of  new  drops  on  the  bare  area  is  continued  for  stage 
two  until  the  termination  criteria  is  met.   At  this  time,  the 
third  begins  with  a  larger  area  and  uses  a  starting  process 
which  is  analogous  to  the  start  up  of  the  second  stage. 
Time  steps  equal  one  half  the  cycle  time  of  the  second  stage 
are  used.   CFigure  2  illustrates  the  start  up  of  the  third 
stage).   During  the  growth  of  the  third  stage,  bare  area  is 
assigned  an  average  heat  transfer  coefficient  found  from  the 
results  of  the  second  stage  for  the  same  time  interval. 
Later  Stages 

Following  stages  are  all  similar  to  the  third 
stage.   They  use  the  same  effective  coefficient  concept  and 
termination  criteria.   Each  stage  differs  from  the  former  by 
having  a  lower  magnification  and  having  a  stage  area  which 
is  ten  times  the  larger  than  the  previous  stage  area.   The 
simulation  finally  ends  on  the  last  stage. 

The  final  stage  area  is  slightly  greater  than  a 
factor  of  ten  larger  than  the  proceeding  stage  area.   The 
size  chosen  is  three  and  one  quarter  square  centimeters  and 
is  comparable  to  the  size  of  the  experimental  condensing 
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surfaces  [2,13].   The  termination  of  the  last  stage  is 
different  from  the  criteria  used  for  the  earlier  stages. 
The  last  stage  shuts  off  when  a  drop  of  the  departing  size 
is  first  formed.   The  values  of  the  departing  drop  size  used 
in  this  simulation  were  experimentally  determined  by  Graham 
[2] .   The  final  output  of  the  complete  simulation  is  a  time 
averaged  heat  transfer  coefficient  for  the  process.   The 
simulation  program  is  visualized  by  a  flow  diagram,  Figure  3. 
For  the  program,  the  input  variables,  source  listing,  and 
description  of  output  are  presented  in  Appendix  C,  Appendix 
D,  and  Appendix  E,  respectively. 
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INPUT 


CALCULATION  OF  CONDENSATION  GROWTH  MODEL 


* — Generation  and  Ordering  of  Sites 
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Increase  Time  by 
Time  Step 
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Test  for  Stage  1 
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no.  of  drops  for  the  two 
step  stage  start  up 


Test  for  Stage  Termination 


No 


Yes 


Commence  a  new  Stage 


No 


Yes 


Stage  Counter  to  test  for 
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FIGURE  3 
PROGRAM  FLOW  DIAGRAM 
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CHAPTER  III 

DISCUSSION  OF  RESULTS 
Effect  of  Site  Density 

At  a  saturation  temperature  of  212 °F  and  at  a 

temperature  difference  of  one-half  degree,  heat  transfer 

predictions  were  made  for  various  site  densities. 

Figure  4  displays  the  time  average  heat  transfer  coefficient 

over  a  complete  simulation  cycle  for  a  range  of  site 

densities.   The  process  is  controlled  by  two  limiting 

factors.   Initially,  the  surface  is  covered  by  small  drops 

which  provide  little  resistance  to  the  heat  flux.   They 

limit  the  process  of  heat  transfer  because  they  are  in 

insufficient  number  to  provide  an  adequate  amount  of  drop 

surface  on  which  vapor  can  condense.   As  time  increases,  the 

average  drop  size  increases  and  the  amount  of  drop  surface 

available  for  condensation  increases.   This  effect  accounts 

for  the  initial  increase  of  the  time  averaged  coefficient. 

The  heat  transfer  coefficient  reaches  a  maximum 
where  another  limiting  factor  becomes  controlling.   As  the 
average  drop  size  increases,  the  average  drop  resistance  to 
heat  transfer  increases.   After  the  maximum,  the  average 
drop  resistance  becomes  the  controlling  limitation  of  the 
process.   At  higher  site  densities,  the  process  is  initially 
less  limited  by  the  inadequate  amount  of  drop  surface  area 
and  becomes  resistance  limited  more  rapidly. 

The  higher  site  densities  enable  the  process  to 
obtain  greater  time  averaged  heat  transfer  coefficients. 


-32- 

Higher  coefficient  values  are  indicative  of  faster  rates  of 
mass  accumulation.   For  this  reason,  drops  of  the  departing 
size  are  produced  sooner  for  process  with  higher  site 
densities. 

Once  the  process  becomes  resistance  limited,  the 
rate  of  decline  in  the  time  average  value  of  the  heat  trans- 
fer coefficient  is  the  same  for  all  site  densities.   This 
observation  is  based  on  the  similarity  of  the  slopes  of  the 
curves  after  the  maximum  has  been  reached.   The  decrease  in 
the  heat  transfer  coefficient  after  the  maximum  appears  to 
be  independent  of  site  density. 

At  the  saturation  temperature  of  88 °F  and  for  the 
temperature  difference  of  one-half  degree,  similar  heat  trans- 
fer predictions  were  made  and  are  plotted  on  Figure  5.   At 
the  lower  saturation  temperature,  the  interfacial  resistance 
to  heat  transfer  is  greater  by  an  order  of  magnitude  and  has 
the  effect  of  decreasing  the  growth  rate  of  smaller  drops  by 
a  similar  amount.   The  reduced  growth  rate  at  lower  satura- 
tion temperatures  makes  the  process  more  strongly  limited 
by  the  initial  inadequacy  of  available  condensing  surface. 
Consequently,  at  the  lower  saturation  temperature  smaller 
maximums  occur  later  in  time  when  the  site  density  is  held 
constant. 

A  calculation  of  the  maximum  theoretical  heat 
transfer  coefficient  based  on  a  maximum  drop  growth  rate 
was  made.   The  details  of  the  calculation  are  given  in 


-33- 


Appendix  F.   The  values  obtained  are  plotted  on  Figure  6 
which  is  a  graph  of  the  time  averaged  heat  transfer 
coefficient  at  the  end  of  the  cycle  for  different  site 
densities.   The  site  densities  examined  produced  coefficients 
well  below  the  theoretical  maximum.   For  a  given  site  den- 
sity, a  larger  value  of  the  coefficient  was  obtained  at  the 
higher  saturation  temperature  because  of  the  faster  growth 
rates.   As  the  site  density  increases,  the  difference  in 
growth  rates  has  a  more  pronounced  effect  and  the  coeffi- 
cient values  for  the  two  saturation  temperatures  diverge. 
As  the  site  density  decreases,  the  effect  of  the  difference 
in  growth  rates  diminishes  and  the  coefficient  values 
converge. 

Next,  the  time  averaged  distribution  of  drop 
sizes  over  the  complete  cycle  for  various  site  densities 
were  plotted  and  compared.   Figure  7  and  Figure  8  show  the 
number  of  drops  per  square  centimeter  for  a  given  drop 
radius  at  the  saturation  temperature  of  88°F  and  212 °F.   At 
either  saturation  temperature  and  for  drop  radii  greater  than 
ten  microns,  the  distributions  are  almost  invariant  with 
site  density.   For  radii  smaller  than  ten  microns,  the 
distributions  diverge.   At  the  higher  site  densities,  there 
are  more  small  drops.   The  presence  of  greater  numbers  of 
small  drops  account  for  the  larger  heat  transfer  coefficients 
predicted  at  higher  site  densities. 
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Estimation  of  Experimental  Site  Densities 

A  determination  of  actual  site  density  was  made  by 
comparirg  the  model's  drop  distribution  for  the  simulation 
cycle  to  Graham's  measured  distribution  of  the  steady  state 
process.   Figures  9  and  10  display  the  model's  prediction 
and  Graham's  data  points  for  approximately  the  same  heat 
transfer  coefficient.   The  theoretical  distributions  agree 
well  with  Graham's  measured  distributions.   Since  the  model's 
distributions  are  invariant  with  site  density  in  the  region 
measured  by  Graham,  the  actual  site  density  had  to  be 
inferred  by  matching  the  model's  predicted  coefficient  with 
Graham's  measured  heat  transfer  coefficient.   For  a  tempera- 
ture difference  of  one-half  degree,  the  model  indicated  that 

Graham's  smooth  copper  condensing  surface  had  an  active  site 

7 
density  of  about  2.5  x  10   sites  per  square  centimeter  at 

212°F  and  of  about  2.0  x  106  at  88°F. 

A  standard  experimental  method  of  determining  site 
density  is  to  photograph  a  portion  of  the  condensing  surface 
with  a  large  number  of  just  visible  drops.   These  drops  have 
diameters  close  to  the  wave  length  of  visible  light;  there- 
fore, the  picture's  quality  is  limited  by  the  resolution  of 
optics.   By  estimating  the  size  of  the  drops  and  their  num- 
ber, a  drop  density  can  be  calculated.   The  nucleation  site 
density  can  be  assumed  to  be  equal  to  or  greater  than  the 
measured  drop  density.   The  model  indicated  that  this 
technique  may  yield  an  incorrect  estimate  at  high  site 
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densities.   For  various  nucleation  site  densities  of  the 
model,  Table  2  shows  the  model's  highest  instantaneous 
drop  density  for  different  radii.   Even  if  drops  with  a 
radius  of  one  micron  could  be  accurately  measured,  the 
nucleation  site  density  could  not  with  certainty  be 
inferred. 


10' 
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FIGURE  10   COMPARISON  OF  PREDICTED  AND  MEASURED 
DROP  DISTRIBUTIONS  AT  212 °F 


-43- 


TABLE  2 
HIGHEST  INSTANTANEOUS  DROP  DENSITY 


Model's  Nucleation  Drop  Radius  in  Microns 

Site  Density        1.0        2.5        5.0       10.0 

106  8.1  x  105  7.4  x  105   3.6  x  105   9.0  x  104 

107  3.6  x  106  1.4  x  106   3.3  x  105   6.0  x  104 

108  7.6  x  106  2.3  x  106   2.2  x  105   9.0  x  104 

109  1.0  x  107  2.0  x  106   3.7  x  105   8.0  x  104 


Note:   The  site  densities  are  in  dimensions  of  sites  per 
square  centimeter.   The  drop  densities  have  units 
of  drops  per  square  centimeter.   The  predictions 
are  for  the  condensing  conditions  of  one-half 
degree  temperature  difference  and  for  a  saturation 
temperature  of  212 °F. 
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Effect  of  Temperature  Differences 

At  the  saturation  temperature  of  88 °F  and  212 °F  and 
for  the  site  densities  determined  by  matching  Graham's 
measured  heat  transfer  coefficients,  predictions  were  made 
for  various  temperature  differences  between  the  vapor  and 
the  surface.   The  results  are  displayed  in  Figure  11  and 
Figure  12.   The  model  predicts  an  approximately  constant 
heat  transfer  coefficient  for  a  constant  site  density  in  the 
temperature  difference  range  from  a  quarter  degree  to  four 
degrees .   The  experimental  results  also  plotted  on  the 
figures  agree  with  the  constant  heat  transfer  coefficients 
values  for  temperature  differences  greater  than  one  degree. 
Below  one  degree  the  heat  transfer  coefficients  are 
experimentally  observed  to  decrease.   This  disagreement 
between  the  model  prediction  for  low  temperature  difference 
may  be  attributable  to  the  model's  assumption  of  constant 
site  density. 

The  minimum  drop  radius  increases  rapidly  as  the 
temperature  difference  is  reduced.   Figure  13  is  a  graph  of 
the  minimum  drop  radius  plotted  as  a  function  of  temperature 
difference.   The  curve  indicates  that  as  the  temperature 
difference  decreases  from  one  degree  to  a  quarter  degree, 
the  minimum  drop  radius  increases  by  a  factor  of  five  for 
both  88 °F  and  212°F  saturation  temperatures.   At  high  site 
densities,  large  values  of  the  minimum  radius  could  cause 
artificial  overlapping  on  the  first  stage  of  the  model  and 
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thereby  predict  an  erroneously  high  heat  transfer  coefficient 

The  maximum  theoretical  site  density  for  a  given 
temperature  difference  was  determined  to  see  if  the  con- 
stant site  density  assumption  at  small  temperature 
differences  was  causing  artificial  overlapping.   Appendix  G 
contains  derivation  of  a  maximum  site  density  for  a  given 
temperature  difference.   Figure  14  is  the  result  of  a  calcula- 
tion using  the  derivation.   Even  for  a  temperature  difference 
of  a  quarter  degree,  the  theoretical  maximum  density  is  well 
above  the  site  densities  assumed  for  this  analysis.   The 
model's  predictions  at  temperature  differences  less  than  one 
degree  are  believed  to  be  valid  for  an  assumption  of  constant 
site  density.   Kinematics  of  drop  growth  cannot  account  for 
the  measured  decrease  in  the  heat  transfer  coefficient. 
The  observed  decrease  can  be  explained  by  a  reduction  in  the 
actual  site  density  as  the  temperature  difference  decreases 
below  one  degree. 
Effect  of  Saturation  Temperature 

To  isolate  the  effect  of  saturation  temperature, 
calculations  were  made  for  the  saturation  temperatures  of 
88°F,  150°F/  212°F,  and  274°F  at  the  site  density  of  107  and 
at  the  temperature  difference  of  three  degrees.   Figure  15 
displays  the  theoretical  prediction  for  the  heat  transfer 
coefficient  as  a  function  of  saturation  temperature  at  the 
constant  site  density.   Above  150 °F  the  coefficient  remains 
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approximately  constant.   Below  150 °F  the  heat  transfer 
coefficient  decreases. 

The  cause  of  the  decrease  can  be  understood  by 
examining  the  drop  distributions  for  the  four  saturation 
temperatures.   Figure  16  displays  the  drop  distributions. 
As  the  saturation  temperature  increases ,  the  number  of 
small  drops  increases.   This  effect  can  be  explained  by 
considering  the  interfacial  resistance.   As  the  saturation 
temperature  increases,  the  interfacial  resistance  decreases; 
and  the  drop  growth  rate  of  the  smaller  drops  increases . 
For  the  constant  site  density  assumption,  it  is  the  increased 
growth  rate  that  is  responsible  for  the  greater  numbers  of 
small  drops  which  appear  in  the  drop  distribution.   Higher 
heat  transfer  coefficients  are  predicted  when  the  drop 
distribution  contains  more  small  drops. 

Graham's  experimental  data  for  heat  transfer 
coefficients  as  a  function  of  saturation  temperature  is  also 
plotted  on  Figure  15.  The  experimental  result  has  a  greater 
slope  than  the  constant  site  density  curve.  This  difference 
of  slopes  is  further  indication  that  the  actual  site  density 
decreases  as  the  saturation  temperature  decreases  from  212 °F 
to  88°F. 
Time  Necessary  to  Produce  a  Drop  of  the  Departing  Size 

At  a  saturation  temperature  of  212 °F  and  for  a 
temperature  difference  of  two  degrees,  predictions  were  made 
by  varying  site  density  in  order  to  match  Westwater's  data  [13] 
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for  heat  transfer  coefficients  and  sweeping  frequency  of 
drops  of  the  maximum  size.   The  frequencies  were  converted 
into  the  time  necessary  to  produce  a  drop  of  the  departing 
size.   Figure  17  shows  measured  and  predicted  heat  transfer 
coefficients  plotted  as  a  function  of  the  time  necessary  to 
produce  departing  drops.   For  the  same  heat  transfer 
coefficient,  the  model  predicts  times  which  are  twice  as 
large  as  Westwater  observed.   Figure  18  shows  the  model's 
time  averaged  coefficient  as  a  function  of  simulation  time 
at  the  same  condensity  conditions.   Reducing  the  model's 
cycle  time  by  a  factor  of  two  means  only  a  twenty  percent 
increase  in  the  model's  predicted  coefficient.   For  this 
reason,  the  lack  of  agreement  with  measured  times  is  not 
considered  to  invalidate  the  other  results  of  the  model. 

This  disagreement  may  be  attributed  to  the 
hemispherical  drop  assumption  or  to  the  model's  failure  to 
account  for  sweeping.   Drops  of  all  sizes  were  assumed  to  be 
hemispherical  in  shape.   A  flatten  or  distorted  drop  of  the 
same  volume  as  a  hemispherical  drop  will  cover  a  larger 
portion  of  the  surface.   The  larger  surface  coverage  for  the 
same  drop  mass  implies  a  greater  opportunity  for  coalescence 
Since  large  drops  grow  mostly  by  coalescence  [2,11],  the 
distorted  drops  should  mature  faster  than  the  perfectly 
hemispherical  drops. 
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The  simulation  is  a  model  of  the  start  of  condensing 
on  an  initially  bare  surface  and  does  not  include  the  effect 
of  sweeping.   Experimentally,  drops  of  the  departing  size 
slide  down  the  surface  and  sweep  vertical  strips  of  the  sur- 
face bare.   In  these  bare  regions,  drops  grow  and  appear 
approximately  uniform  in  size.   Since  a  drop  can  depart  at 
any  time  in  the  actual  steady  state  process,  adjacent  strips 
contain  uniform  drops  of  different  generations.   The  regions 
of  drops  of  different  generations  interact  and  may  hasten  the 
production  of  drops  of  the  departing  size.   The  model's 
failure  to  account  for  sweeping  or  the  hemispherical 
assumption  may  be  the  reason  for  the  overestimate  of  the 
time  necessary  to  produce  a  drop  of  the  departing  size. 
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CHAPTER  IV 
CONCLUSIONS 

1.  Graham's  model  of  heat  flux  through  a  single  drop 
is  a  valid  representation  of  the  physical  phenomenon. 

2.  As  the  site  density  increases,  the  heat  transfer 
coefficient  increases  because  of  a  greater  number  of  small 
drops  in  the  distribution. 

3.  For  a  temperature  difference  greater  than  one 
degree,  the  measured  heat  transfer  coefficient  remains 
constant  because  the  site  density  is  constant.   For  tempera- 
ture differences  less  than  one  degree,  the  decrease  in 
measured  coefficients  can  be  explained  by  a  reduction  in  the 
site  density. 

4.  As  the  saturation  temperature  decreases  from  212°F 
to  88°F,  the  active  nucleation  site  density  decreases. 

5.  The  model  predicts  a  time  necessary  to  produce  a 
drop  of  the  departing  size  which  is  twice  as  large  as  is 
observed  for  the  steady  state  process. 

6.  This  computer  simulation  provides  a  method  of 
accurately  predicting  the  heat  transfer  properties  of  the 
dropwise  condensation  process. 
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CHAPTER  V 
SUGGESTIONS  FOR  FUTURE  WORK 

1.  The  present  model  can  be  improved  by  using  finer 
time  steps  and  more  elaborate  starting  technique  for  the 
later  stages.   Incorporation  of  the  sweeping  effect  should 
also  be  added  so  that  the  simulation  will  more  closely  model 
the  steady  state  process. 

2.  The  effect  of  non-hemispherical  drops  may  require 
investigation.  A  combination  of  experimental  observations 
and  a  theoretical  analysis  should  determine  the  effect, 

if  any,  of  non-hemispherical  drops. 
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APPENDIX  A 
MODEL  OF  DROP  GROWTH  DUE  TO  CONDENSATION 

Graham  [2]  proposed  an  analytical  model  for  the 
growth  of  a  single  drop  due  to  condensation  of  vapor.   The 
heat  flux  associated  with  the  condensing  of  vapor  passes 
through  three  independent  series  resistances.   The  heat 
transfer  resistances  included  are  curvature,  interfacial, 
and  conduction.   Each  of  these  resistances  is  discussed 
below. 

1.   Curvature  Resistance.   At  a  curved  liquid-vapor 
interface,  the  saturation  temperature  of  the  vapor  in 
equilibrium  with  the  drop  is  less  than  the  saturation 
temperature  of  the  vapor  in  equilibrium  with  liquid  at  a 
flat  interface  at  the  same  pressure.   A  derivation  of  the 
temperature  drop  through  the  interface  due  to  curvature  can 
be  found  in  Graham's  thesis.   The  result  of  his  work  is 
reproduced  below: 

2Ts°   1 

ATc=  h^  I  t*'1' 

where  AT  temperature  drop  due  to  curvature 
T  saturation  temperature 
a  surface  tension 
H~  latent  heat  of  transformation 
p  liquid  density 
r  drop  radius 
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When  the  total  temperature  difference  available  between  the 
vapor  and  the  condensing  surface  is  equal  to  the  temperature 
drop  due  to  curvature,  the  minimum  drop  size  which  can 
exist  on  the  surface  is  determined.   The  expression  for 
minimum  drop  radius  is  written  below: 

2V     1 
rmin  "   Hf  p   ATfc  (A,2) 

where  AT.  the  total  temperature  difference  between 
the  vapor  and  the  condensing  surface. 
By  combining  the  two  previous  equations ,  an  expression  for 
the  temperature  difference  due  to  curvature  can  be  found. 

ATc  =   -2^2:  ATt  (A.  3) 

where  AT  the  temperature  difference  due  to 
curvature . 
This  temperature  difference  is  not  an  actual  resistance  to 
heat  transfer,  but  can  be  thought  of  as  an  equivalent 
resistance. 

2.   Interfacial  Mass  Transfer  Resistance.   At  the  liquid- 
vapor  interface,  a  pressure  difference  is  necessary  to  force 
the  mass  transfer.   The  pressure  difference  can  be  converted 
into  a  temperature  difference  and  is  written  below: 


AT.  =  — S — j  (A.  4) 

hi2irr 
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Where  Q  rate  of  heat  transfer 

h.  interfacial  heat  transfer  coefficient 

AT.  temperature  drop  due  interfacial 
resistance. 

The  interfacial  coefficient,  h . ,  can  be  calculated  from  the 

equation  below  which  was  derived  by  Nabavian  and  Bromley  [14] . 

H  -    &S>     ^->V2  A-  CA.3) 

X  *   a    2TTRT         s  Vg 

s 

where  a  condensation  coefficient  (assumed  to  be 
unity) . 

M  molecular  weight 
R  universal  gas  constant 
V  specific  volume  of  the  vapor. 
Figure  A-l  is  a  plot  of  the  interfacial  heat  transfer  coeffi- 
cient as  a  function  of  saturation  temperature. 

3.   Conduction  Resistance.   The  most  significant 
resistance  associated  with  the  liquid  drop  is  the  resistance 
to  heat  conduction  through  the  drop.   The  conduction  resistance 
of 


■o  _   AT  _    r  ,_  -. 

Rcyl    ~Q     r~"2  (A'6) 
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-5 
this  shift.   The  age  discrepancy  is  initially  about  4  x  10 

seconds  and  decreases  as  the  drop  radius  increases .   The 

error  induced  by  the  imprecise  knowledge  of  the  nucleation 

size  should  not  be  significant  because  little  vapor  is 

condensed  in  the  small  time  interval. 

Each  resistance  has  been  expressed  as  a  temperature 
drop.   The  total  temperature  difference  between  the  vapor  and 
the  condensing  surface  is  equal  to  the  sum  of  the  three 
resistance  induced  temperature  drops.   Then  by  utilizing  the 
idea  that  the  rate  of  heat  transfer  through  a  single  drop 
is  proportional  to  the  drop  volume  increase  with  time , 
Graham  obtained  an  expression  for  drop  age  as  a  function  of 
drop  diameter. 

Graham's  derivation  is  reproduced  below  starting 
with  the  summation  of  temperature  drops.  AT.  is  the  total 
temperature  difference  between  the  vapor  and  the  surface. 


ATt  = 


ATC 
(curvature) 

+ 

ATi 
Cinterfacial) 

+ 

ATdc 
(conduction) 

(A. 8) 

min      AT      + 

Q 

4K7rr 

(A. 9) 

r             c          hi2irr 

ATt  = 


The  heat  flux  through  the  drop  is  assumed  to  be  proportional 
to  the  change  in  volume  of  the  drop. 

■  -  p  %  £ 
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where  Q  rate  of  heat  transfer 

AT  temperature  across  the  cylinder 
R   ,  conduction  resistance  of  the  cylinder 
K  thermal  conductivity  of  liquid 
r  drop  radius 
For  hemispherical  drops,  the  above  expression  is  multiplied 
by  a  shape  factor  of  less  than  one.   Mikic  [15]  determined 
the  value  of  the  shape  factor  to  be  one  fourth. 

Multiplying  by  the  shape  factor  and  then  solving 
for  the  temperature  drop  through  the  hemispherical  drop,  an 
expression  for  conduction  resistance  is  obtained: 


AT,   =    Qr,  (A.  7) 


where  AT,  temperature  drop  through  the  drop  due 
to  conduction  resistance. 

Figures  A-2  and  A-3  are  reproduced  from  Graham's  thesis  and 

indicated  the  effect  of  the  three  resistances  on  drop  diameter 

as  a  function  of  drop  age. 

Graham  [2] ,  stated  that  nucleation  radius  was 

probably  between  one  and  two  times  the  minimum  drop  radius. 

A  drop  of  the  minimum  size  has  a  theoretical  zero  growth 

rate.   Several  nucleation  radii  were  examined  to  determine 

their  effect.   The  curve  of  drop  radius  as  a  function  of  drop 

age  kept  essentially  the  same  shape,  but  was  shifted  to  the 

right  as  nucleation  radius  decreased.   Figure  A-4  illustrates 
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where  V  drop  volume 

t  time  or  drop  age 
dV/dt  change  in  drop  volume  with  respect  to  time 


For  an  hemispherical  drop, 


dV  _  .2  dr 

ar  -  27Tr     at 


(A. 10) 


Q  =  pHfg  2™2   H 


(A.  11) 


Substituting  equation  CA.ll)  into  equation  (A. 9), 


at  -  J2HL  at  +  tha.    ^r     ^fg  rdr 
Q1t     r   Ait     h.    dt      2K    dt 


(A. 12) 


Solving  for  dr/dt  and  substituting  drop  diameter  D  for  twice 
the  drop  radius , 


dD  =  -±-     at 

at  ■  pHfg   ait 


1  -  D  .  /D 
mi^ 


D/2K  +  2/h. 


(A. 13) 


Equation  (A. 13)  can  be  integrated  and  then  evaluated  at  the 
lower  limit.   Since  the  logarithmic  terms  become  undefined  at 
drop  age  zero  when  the  drop  diameter  equals  the  minimum 
diameter,  the  function  must  be  evaluated  at  the  minimum 
diameter  plus  epsilon.   Graham's  final  result  of  drop  age  as 
a  function  of  drop  diameter  is  listed  below. 
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~  [".  SCD-Dm.  }2  +  2D.  CD-D  .  }  +  D^    In  CD-D.*] 
2K        mm       min    min     min        min  J 


E7  [CD-Dmin>  +  Dmin  ln  ®-D: 

1     L 


m 


5F  D!^n  "  r  Dmi„  =   tit-  W*-*  (A.  14) 

2K   mm   h.    mm   pH£_    t 


The  simulation  growth  routine  required  an  expression 
of  drop  radius  as  a  function  of  age.   For  all  cases  in  this 
study,  Graham's  equation  (A. 14)  was  inverted  by  substituting 
a  constant  equal  to  five  minimum  drop  diameter  for  the  dia- 
meter variable  in  the  logarithmic  terms.   This  action  con- 
verts the  logarithmic  expressions  into  constants.   The 
inverstion  is  then  finished  by  completing  the  square.   The 
final  result  (A. 15)  is  an  analytical  approximation  for  drop 
diameter  as  a  function  of  drop  age  and  is  shown  below: 


where  a  =  1/4K 

b  "  °min/2K  +  2/h< 

mm'  '  i 

c  =  5  Dmin/4K  '  4Dmir/hi 


d  -  4ATt/pHfg 
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The  approximation  has  reasonable  accuracy  over  almost  all  of 
the  range  of  drop  sizes.   The  greatest  inaccuracy  experienced 
for  the  condensing  conditions  of  this  study  was  about  two 
percent  for  drop  radii  of  less  than  one  micron.   As  the  drop 
radius  increased  above  one  micron,  the  inaccuracy  rapidly 
diminished. 

Using  equations  (A. 14)  and  (A. 15),  the  drop  growth 
routine  can  be  carried  out.   First  a  drop's  age  is  determined 
using  equation  (A. 14).   The  age  calculated  is  then  increased 
by  the  time  step.   Finally,  the  new  size  of  a  drop  is 
calculated  with  expression  CA.15) . 
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APPENDIX  B 

COALESCENCE  ROUTINE 

The  coalescence  routine  developed  is  capable  of 
determining  coalescence  of  drops,  blanketing  of  drop  sites 
or  locations,  and  the  uncovering  of  sites.   Four  variables 
associated  with  each  site  are  required  for  the  routine  to 
operate  with  either  randomly  or  uniformly  distributed  sites. 
The  first  two  of  the  four  site  variables  are  coordinates, 
"X"  and  "Y" .   These  are  obtained  by  superimposing  a  first 
quadrant  coordinate  system  on  the  stage  area.   The  coordin- 
ates of  sites  are  necessary  to  calculate  the  distances  be- 
tween sites.   Then  the  sites  are  ordered  from  left  to  right 
on  the  stage  square  and  assigned  index  numbers  represented 
by  the  variable  "I".   The  site  with  the  lowest  X  coordinate 
value  is  site  number  one;  the  site  with  the  highest  X 
coordinate  value  has  the  largest  index  number.   The  fourth 
variable,  "R",  identifies  the  activity  at  the  site.   The 
"R"  variable  can  assume  three  types  of  possible  values; 
minus  one,  zero,  or  a  finite  positive  number.   The  minus 
one  value  indicates  that  the  site  is  covered  or  blanketed  by 
a  drop  centered  on  another  site.   The  zero  value  represents 
a  site  that  is  exposed  to  the  vapor  or  uncovered.   The  third 
"R"  value  is  a  finite  positive  number  which  is  the  radius  of 
a  drop  in  feet  centered  on  the  site.   Touching  of  drops  and 
covering  of  nucleation  sites  is  determined  by  comparison  of 
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the  distance  between  sites  and  the  drop  radius.   These 
four  site  variables;  "X,  Y,  I,  and  R" ;  are  necessary  to 
perform  the  bookkeeping  required  by  the  coalescence  activity. 

The  coalescence  routine  commences  for  each  time 
step  at  the  site  with  the  lowest  index  number  with  a  drop 
centered  on  it  and  proceeds  until  the  site  with  a  drop  and 
the  highest  index  number  is  examined  for  possible 
coalescences.   A  search  area  around  the  "Ith"  drop  of 
interest  is  defined  as  a  square  of  two  "Ith"  drop  diameters 
on  a  side.   Figure  B-l  is  a  diagram  of  a  search  square  on  the 
simulated  condensing  surface.   Sites  with  a  zero  "R"  value  or 
a  positive  "R"  value  inside  the  square  are  examined.   A 
determination  is  made  whether  or  not  the  "Ith"  drop  with  its 
increased  radius  from  condensation  growth  blankets  any  un- 
covered or  active  sites.   At  the  same  time,  a  check  for  the 
possible  coalescence  between  the  central  "Ith"  drop  and  a 
neighboring  "Jth"  drop  with  a  center  inside  the  square  is 
conducted.   Only  neighboring  drops  of  the  same  size  or 
smaller  are  not  examined  even  though  they  may  be  blanketing 
part  of  the  search  or  even  touching  the  central  drop. 

By  maintaining  a  consistant  examining  procedure  as 
the  routine  moves  from  drop  to  drop  and  from  left  to  right  on 
the  simulated  surface,  all  possible  coalescences  are  found. 
This  approach  requires  a  relatively  small  amount  of  computer 
core  storage  and  is  reasonably  efficient  for  one  thousand 
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sites  on  a  stage  area.   For  large  numbers  of  small  drops, 
small  search  areas  are  used.   For  the  small  number  of  large 
drops,  .large  search  areas  are  employed.   This  flexibility 
is  the  principle  reason  for  the  routine's  efficiency. 

When  a  site  is  found  to  be  covered,  its  "R"  value 
is  changed  from  zero  to  minus  one.   If  the  "Ith"  drop  is 
found  to  be  touching  a  neighboring  "Jth"  drop  inside  the 
search  square,  a  coalescence  occurs.   The  masses  of  the 
"Jth"  and  "Ith"  drops  are  combined  to  form  a  new  drop  at  the 
central  location  of  the  "Ith"  drop.   A  check  for  uncovering 
of  sites,  previously  blanketed  by  the  "Jth"  drop,  is  then 
carried  out  in  a  similar  manner  as  the  coalescence  was 
determined.   When  all  sites  previously  covered  by  the  "Jth" 
drop  have  been  found  and  uncovered,  the  search  for 
coalescence  starts  once  again  by  defining  a  new  larger 
search  square  around  the  "Ith"  drop.   The  increased  size  of 
the  square  is  indicative  the  the  increase  in  size  of  the 
"Ith"  drop  due  to  the  coalescence  with  the  "Jth"  drop. 
When  no  more  coalescences  are  obtained  between  the  "Ith" 
drop  and  any  of  its  neighbors  inside  the  search  square,  the 
coalescence  routine  centers  its  attention  on  the  drop  of 
the  next  higher  index  number.   The  same  search  for 
coalescence,  covering,  and  uncovering  is  carried  out. 
When  the  drop  with  the  highest  index  number  fails  to  touch  a 
neighboring  drop,  the  coalescence  routine  ends. 
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APPENDIX  C 
DESCRIPTION  OF  THE  INPUT  DATA 

Seven  data  cards  are  used  to  read  in  fifteen 
variables  necessary  to  run  the  program.   These  variables  will 
be  discussed  in  the  order  which  they  are  read  in. 

1.  First  Data  Card 

Variable  Name :   IX 

Format:   19 

Description:   The  variable  "IX"  is  any  number  with 
nine  or  less  integers  so  long  as  the  last  integer  is  odd. 
The  number  is  necessary  to  call  a  random  number  generator 
which  is  the  IBM  Scientific  Subroutine  RANDU.   The  random 
number  generator  is  used  to  create  the  random  site  location 
for  the  first  stage.   If  fewer  than  one  hundred  sites  are 
used  for  the  first  stage,  the  selection  of  the  variable  "IX" 
can  cause  the  time  averaged  heat  transfer  coefficient  at  the 
end  of  the  first  stage  to  vary  by  as  much  as  ten  percent. 

2.  Second  Data  Card 

Variable  Names:   NN,  SCALE 

Format:   110,  F10.4 

Description:   These  two  variables  determine  the  site 
density.   "NN"  is  the  number  of  sites  for  the  first  stage. 
"SCALE"  is  the  dimension  of  the  side  of  the  first  stage 
square.   Scale  is  read  in  units  of  microns  and  then  conver- 
ted into  feet  in  the  program.   (The  program  does  most  of  its 
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calculations  in  English  Engineering  Units.)   Any  number  of 
"NN"  sites  can  be  used,  but  the  scale  can  be  only  10.0, 
31.623,  100.0,  316.23,  or  1000.0  microns.   Site  densities  up 
to  ten  to  the  ninth  sites  per  square  centimeter  can  be 
accommodated.   (The  program  ran  more  efficiently  by  not 
letting  the  number  of  sites  exceed  one  thousand  on  the 
initial  stage  area.) 

3.  Third  Data  Card 

Variable  Names:   MM,  TSTEP 

Format:   110,  Ell. 4 

Description:   "MM"  is  the  number  of  time  steps 
allowed  for  each  stage.   A  limit  of  one  thousand  should  not 
be  exceeded.   "TSTEP"  is  the  time  step  for  the  first  stage. 
The  value  selected  should  be  large  enough  so  that  drops 
originating  from  nucleation  sites  have  a  radius  which  is 
three  or  more  times  larger  than  the  minimum  drop  radius.   The 
growth  models  inaccuracy  for  very  small  drops  dictates  this 
constraint.   TSTEP  has  the  units  of  seconds. 

4.  Fourth  Data  Card 

Variable  Names:   TSAT,  SURTE,  HEATL,  THCON,  DTEMT, 

FDEN,  HTCI 
Format:   6F10.4,  Ell. 4 
Description: 

"TSAT"  is  the  saturation  temperature  of  the 
vapor  for  a  flat  vapor-liquid  interface  in  degrees  Rankine. 
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"SURTE"  is  the  surface  tension  of  the  liquid 
at  the  saturation  temperature  and  has  the  units  of  pounds 
force  per  foot. 

"HEATL"  is  the  latent  heat  of  transformation 
of  the  vapor  to  the  liquid  at  the  saturation  temperature. 
Its  units  are  British  Thermal  Units  per  pound  mass. 

"THCON"  is  the  thermal  conductivity  of  the 
liquid  at  the  saturation  temperature  and  has  the  units  of 
British  Thermal  Units  per  foot-hour-degree  absolute. 

"DTEMT"  is  the  temperature  difference  between 
the  vapor  and  the  base  of  the  drop  which  is  approximated  by 
the  average  surface  temperature.   The  units  are  degrees 
Rankine . 

"FDEN"  is  the  liquid's  density  in  pounds  mass 
per  cubic  foot. 

"HTCI"  is  the  inter facial  heat  transfer 
coefficient  and  was  obtained  from  a  graph  of  values  for 
various  saturation  temperature  in  Graham's  thesis  [2] .   The 
coefficient  has  the  units  of  British  Thermal  Units  per 
square  foot-hour-degree  absolute. 
5.   Fifth  Data  Card 

Variable  Name:   RNUC 
Format:   F6.3 

Description:   "RNUC"  is  the  size  chosen  for  the 
nucleating  drops  necessary  for  the  math  model  of  condensation 
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growth.   Graham  [2]  stated  that  the  nucleation  size  was 
probably  between  one  and  two  times  the  minimum  diameter. 
This  simulation  used  a  nucleation  drop  radius  of  one  and  a 
half  times  the  minimum  radius.   The  "RNUC"  value  for  this 
case  is  the  factor  one  point  five. 
6.   Sixth  Data  Card 

Variable  Name:   RMAX 

Format :  F10 . 4 

Description:   "RMAX"  is  the  radius  in  microns  of  a 
drop  of  the  departing  size. 

The  values  used  in  this  study  were  determined  by 
Graham  [2] .   For  a  temperature  difference  of  one-half 
degree,  Graham  measured  values  of  1500  microns  and  1250 
microns  for  the  saturation  temperatures  of  88  and  212 
degrees  Fahrenheit  respectively.   Values  of  1125  and  1375 
microns  were  assumed  for  the  saturation  temperatures  of 
274  degrees  and  150  degrees.   (For  this  study,  the  maximum 
drop  size  was  assumed  to  remain  constant  for  all  temperature 
differences . ) 
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APPENDIX  D 
PROGRAM  SOURCE  LISTING 

The  computer  simulation  of  this  study  consists 
of  a  single  program  written  in  Fortran  IV,  Level  G,  for 
use  on  the  IBM  Model  360/65  Digital  Computer.   The  follow- 
ing pages  of  this  appendix  contain  a  source  listing  of  the 
program. 
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APPENDIX  E 
DESCRIPTION  OF  OUTPUT 

This  appendix  is  added  to  assist  a  possible  future 
user  of  the  program.   The  output  will  be  categorized  and 
discussed  in  the  order  which  it  appears. 

1.  The  input  variables  defined  in  Appendix  C  are 
written  as  the  first  output.   These  units  are  given.   Next, 
the  result  of  a  calculation  of  minimum  drop  radius  and  the 
variables  for  the  drop  growth  model  are  written  out.   This 
information  is  contained  in  the  first  twelve  lines  of 
output . 

2.  Each  time  step  has  ten  lines  of  output.   The  time 
step  output  variables  are  written  below  as  they  appear. 

"HTC"  is  the  instantaneous  heat  transfer 

2 

coefficient  for  the  time  step  in  BTU/HR  FT  °F. 

"TIME"  is  the  simulation  clock  time  in  seconds. 

"THTC"  is  the  time  averaged  heat  transfer 
coefficient. 

"COALESCENCE  COUNTER"  is  the  number  of  coalescences 
which  occurred  during  the  time  step. 

Distribution  categories  are  listed.   "COVERED" 
implies  that  the  drop  site  is  blanketed  or  covered. 
"PRIMARY"  implies  that  a  drop  site  is  uncovered.   Other 
distribution  categories  are  listed  as  drop  radii  in  microns. 
A  radius  category  is  defined  by  a  radius  differential  of 
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forty  percent.   Beneath  the  distribution  categories,  the 
numbers  of  each  kind  is  listed. 

"TOTAL  NO.  OF  DROPS"  is  simply  the  total  number  of 
drops  on  the  stage  surface  at  the  end  of  the  coalescence 
routine. 

"AVERAGE  DROP  SIZE"  is  the  average  size  of  a  drop 
existing  on  the  surface  at  the  end  of  the  coalescence  routine 

"DAREAM"  is  the  amount  of  area  covered  by  drops 
after  coalescence  and  has  the  units  of  square  microns . 

"PERC"  is  the  percentage  of  the  stage  area  covered 
by  drops  of  the  category  which  can  cause  termination  of  the 
stage. 

"BAREAM"  is  the  amount  of  bare  area  in  square 
microns  after  coalescence.   This  variable  is  the  last  of  the 
time  step  variables. 

3.  The  third  major  kind  of  output  data  described  is 
the  start  up  procedure  between  the  stages.  The  variables 
will  be  listed  and  described  as  they  appear. 

"SCALE"  is  the  dimension  of  the  side  of  the  next 
stage  square  in  microns. 

"TSTEP"  is  the  time  step  for  the  next  stage  in 
seconds. 

"THTC1  and  THTC2"  are  the  time  averaged  heat 
transfer  coefficients  obtained  from  the  pervious  stage 
result.   The  start  up  iteration  matches  these  values  for  the 
next  stage. 
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Iteration  sequence  usually  consist  of  several  lines 
on  which  the  iteration  variables  repeated.  "LC"  is  the  number 
of  the  iteration  trial.   "NR1"  is  the  number  of  drops  to  be 
located  on  the  stage  area  for  the  first  time  step.   "RMM1"  is 
the  radius  of  uniform  drops  produced  on  the  bare  area.   "HCK" 
is  a  comparison  of  the  time  averaged  heat  transfer  coeffi- 
cients for  the  two  consecutive  stages  at  the  second  time  step 
of  the  later  stage.   "RMM12"  is  the  radius  in  microns  which 
the  uniform  drops  obtain  after  growing  for  one  time  step. 
The  iteration  ends  when  the  MHCK"  value  is  within  acceptable 
tolerances.   The  number  of  sites  on  the  uniformly  spaced 
stage  area  is  indicated  by  the  variable  "NNLT" .   The 
appearance  of  this  variable  implies  the  completion  of  the 
start  up  process. 

4.  At  the  end  of  the  last  stage,  the  output  for  the 
last  time  step  contains  the  time  averaged  heat  transfer 
coefficient  for  the  entire  process  (THTC) .   The  time  to 
produce  a  drop  of  the  departing  size  is  listed  as  the 
variable  "TIME".   Also  the  size  and  number  of  departing  drops 
are  listed. 

5.  The  last  significant  output  is  an  array  containing 
the  information  of  the  drop  distribution.   The  output 
variables  are  described  as  they  appear. 

"WSSCD"  values  are  the  radius  categories  of  the 
drop  distribution. 
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"BWSS(I)"  values  are  the  time  average  distribution 
of  drops  resulting  from  the  separate  stage  calculations. 
The  units  are  drops  per  square  centimeter. 

"DELN(I)"  values  are  the  time  average  drop 
distribution  for  the  complete  simulation  cycle  and  have  units 
of  number  of  drops  per  square  centimeters . 

"HTCD(I)"  is  the  cumulative  heat  transfer 
coefficient  for  all  drops  of  a  given  size  ("WSS(I)")  or 
smaller. 

"THTCD"  is  the  time  averaged  heat  transfer 
coefficient  for  the  entire  cycle  and  is  calculated  from  the 
distribution  of  drops.   This  output  variable  is  the  last  of 
the  program. 
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APPENDIX  F 
MAXIMUM  THEORETICAL  HEAT  TRANSFER  COEFFICIENT 

An  estimate  of  the  maximum  theoretical  heat 
transfer  coefficient  can  be  made  from  the  maximum  rate  of 
heat  transfer  through  a  single  drop.   Starting  with  an 
expression  for  the  heat  flux  through  a  single,  drop,  the 
maximum  theoretical  coefficient  is  derived  below. 


Q  =  PHfg  g  (F.l) 


For  a  hemispherical  drop, 


dV   -    7-nr2     dr  (V     0\ 

ar   -    27rr    ar  CF-2) 


The  heat  transfer  coefficient  is  defined  as  follows 


h  _  Q   1 

n  "  A  If  (F.3) 

Substituting  equation  (F.3)  into  equation  (F.4), 

h  =  Sl  2«2  If  (F.4, 

2 
The  area  of  the  base  of  the  drop  is  irr  . 

.    pHfg   0dr  ,_  -, 

Substituting  in  dD/dt  for  2dr/dt, 

h  -  PHf^   dD  IF    G) 

h  -  tt   ar  CF*6) 
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The  expression  (F.6)  for  the  heat  transfer  coefficient 
becomes  a  maximum  when  the  drop  growth  rate  variable 
(dD/dt)  obtains  its  maximum  value.   Drop  growth  rate  is 
defined  by  equation  (A. 13)  in  Appendix  A.   Finally/  the 
maximum  theoretical  heat  transfer  coefficient  is  expressed 
below  as  equation  (F.7). 

.       pHfg   ,dD,  ,_  _* 

hmax  =  IT   Wmax  (F'7) 
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APPENDIX  G 
DERIVATION  OF  MAXIMUM  NUCLEATION  SITE  DENSITY  * 


To  ensure  that  the  site  densities  used  in  this 
study  were  realistic,  a  derivation  of  a  theoretical  maximum 
density  was  made.   It  should  be  noted  that  the  theoretical 
maximum  is  much  greater  than  could  be  attained  on  a  real 
condensity  surface. 

Nucleation  centers  are  assumed  to  be  arranged  in 
a  triangular  array,  and  the  distance  between  centers  is 
equal  to  twice  the  minimum  drop  radius  (see  Appendix  A) . 
For  a  triangular  array,  ninety-one  percent  of  the  condensing 
surface  is  covered  if  the  sites  are  circular.   An  equation 
can  be  derived  which  gives  the  maximum  density  as  a  function 
of  the  temperature  difference  and  other  properties. 

918  2 

Maximum  Site  Density  =  - — j   x  10   (sites/cm  ) 

Trr 

CG.l) 

The  radius,  r,  is  equal  to  the  minimum  drop  radius  in 
microns  and  is  defined  below. 

_  2V     i 

rmin   Hf  p   AT  (G.2) 

Substituting  in  equation  (G.2)  into  equation  (G.l), 
the  maximum  theoretical  site  density  as  a  function  of  the 
temperature  difference  and  other  properties  can  be  obtained. 
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